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Two-photon decay rates of pseudoscalar, scalar and tensor states of charmonium and 
bottomonium are calculated in the framework of the relativistic quark model. Both 
relativistic effects and one-loop radiative corrections are taken into account. The obtained 
results are compared with other theoretical predictions and available experimental data. 
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The investigation of the two-photon decay rates of heavy quarkonia is important 
for understanding the heavy quark dynamics in mesons. These decays are sensitive 
probes of the quarkonium wave function. In the nonrelativistic limit their rates 
are proportional to the square of the wave function or its derivative at the origin. 
However, relativistic effects turn out to be substantial, especially for charmonium, 
and significantly modify this dependence. In this note we consider the two-photon 
decays in the framework of the relativistic quark model. 
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Fig. f . Two-photon annihilation diagrams of the quarkonium. 



The diagrams of the QQ bound state with mass M (Q = c, b) annihilation into 
two photons with momenta fci, &2 and polarization vectors e±, £2 are graphically 
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presented in Fig. 1. The corresponding amplitude in the rest frame of the decaying 
meson is then given by 

d 3 p 



7W 77 = 4W3e| yj4L^«£(-p) 
+ (1 «-» 



(p - fci) 2 - m 2 Q 
<,(pin'(P)- (1) 



where a = e 2 /(47r) is the fine structure constant; e<g, mg and p are the charge (in 
units of the electron charge e), mass and relative momentum of quarks; y/3 is the 
colour factor and the Dirac spinors of quark and antiquark arc defined by 
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with p° = eg (p) = ^/p 2 + mg. The photon momenta in the chosen frame obey the 

relations: ki = — k 2 and = k® = |ki| = |k 2 | = M/2. The meson wave function 
for the state with the orbital momentum L, spin S and total momentum J has the 
form 

*m(p) =^(i5mA|JMj) (^st\SX)Y Lm (0, l p)4 >M (p), (3) 

m, A 

where 4>m(p) = </>m(|p|) is the radial wave function. 

For the calculation of the heavy quarkonium wave functions we use the rela- 
tivistic quark model based on the quasipotential approach in quantum field theory. 
In this modelP^lsil a meson is described by the wave function of the bound quark- 
antiquark state, which satisfies the quasipotential equation of the Schrodinger type 
in the center-of-mass frame: 

■b^M) _ jT\ = I A v(p>q;M) ^ (q)> (4) 

2p,R ZfiRj J {2tt)° 

where the relativistic reduced mass is 

M 4 _ (m 2 _ w 2)2 

m = 4M^ > W 

and b (M) denotes the on- mass-shell relative momentum squared 
,2 , m= [M 2 - (mi + m 2 ) 2 ][Af 2 - (m, - m 2 ) 2 ] 

W 4M 2 ' 

Here mi i2 and M are quark masses and a meson mass, respectively. For the QQ 
bound system (quarkonium) mi = m 2 = wlq and Eqs. J^J, © take the form 

M* = Y> b 2 (M) = ^-m* Q . (7) 
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The kernel V(p,q;M) in Eq. JIJ is the quasipotential operator of the quark- 
antiquark interaction. It is constructed with the help of the off-mass-shell scatter- 
ing amplitude, projected onto the positive energy states. The quasipotential for the 
heavy quarkonia, including retardation and one-loop radiative corrections, and ex- 
panded in p^/tHq, can be found in Refs. EEI As the initial approximation in the 
configuration space we obtain the Cornell potential: 



V( r ) = --— + Ar 
3 r 



B, 



(8) 



where ay is the effective coupling constant: ay = 0.450 for charmonium and ay = 
0.285 for bottomonium, respectively. In our model the quark masses = 4.88 GeV, 
m c = 1.55 GeV, m s — 0.50 GeV, m U: d = 0.33 GeV, the slope of the linear potential 
A = 0.18 GeV 2 and the constant term B = —0.16 GeV have fixed values which 
agree with the usually accepted ones. In the recent papeJ^ we calculated the mass 
spectra of heavy quarkonia with complete account of the relativistic corrections of 
order p 2 /mg, including retardation effects, and one-loop radiative corrections. The 
obtained results agree with experimental data within few MeV. Now we apply the 
found relativistic wave functions to the calculation of the two-photon decay rates 
of heavy quarkonia. 

The two-photon decay amplitudes Q of the heavy quarkonium pseudoscalar 
(^o), scalar ( 3 Po) and tensor ( 3 P2) states have the following structure 

Myy^So -> 77) = gtuvpcrk^kZefcZ, 

M J7 ( 3 P -> 77) = 5o[(ei • £a)](fci • ks) - (ei ■ k 2 )(e 2 ■ h)], 
M 17 ( 3 P 2 -» 77) = giKh ■ k^el - (ex ■ k 2 )k{e v 2 

~{^-k 1 )k^e'i + {e 1 -e 2 )k^ 2 ]^ 

+92[(£i • e 2 )(*i ■ k 2 ) - ( £l • k 2 )(e 2 ■ h)]^^" /(h ■ fc 2 ), (9) 

where <?, are corresponding decay constants and is the polarization tensor of 
the tensor meson. Performing angular integrations, using the angular dependence 
of the wave functions given by Eq. (j3J , and summing over photon polarizations 8, we 
get the following relativistic expressions for the decay ratefl 
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d 3 p 



(2tt) 3 Mp 



P m £ Q(p)+P 



(?) cq(p)-P 



M 



^ (2 fgM i n e Q(p) +p 
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>(p) 



(11) 



a For decays of tensor 3 P2 states it is also necessary to average over states with different projection 
Mj of the total momentum J. 
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Similar relativistic expressions have also been found in Refs. |SJ 

In the nonrclativistic limit p/ttlq — > and M — ► 2mg these expressions reduce 
to the known relations^ 



r NR (^ - 77) 

r NR ( 3 Fo - 77) 

r NR ( 3 P 2 -> 77) 



rn 

„2„4 



27a 2 e 4 Q |i?; p (0)r 



36aV Q K P (0)r 

5m Q 



(14) 
(15) 



The first order QCD radiative corrections to the two-photon decay rates can be 
accounted for by multiplying the relativistic decay rates T R pi)|l - (|12(l by the corre- 
sponding factorJS 



(16) 
(17) 
(18) 
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For the numerical analysis we use the two-loop expression for the QCD coupling 
constant a s with = 175 MeV obtained^ from the experimental ratio of T(T — > 
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ggj)/T(T — > ggg). This yields a s (rrib) = 0.18 for bottomonium and a s (m c ) = 0.26 
for charmonium, respectively. 

Table 1. Two-photon decay rates of heavy quarkonium pseudoscalar ( So), scalar ( 3 Po) an d tensor 
( 3 P 2 ) states (in keV). Comparison of different theoretical predictions and experimental data. 



Particle Theory Experiment 





our 


RefGD 


Ref0 


ReflTQl 


RefHTl 


RefH2l 


PDG JS1 


E83^H1 


BellJGD 




5.5 


3.5 


10.94 


7.8 


5.5 


4.8 


7.5(8) 


4.4(1. 7)(1. 5) 






1.8 


1.38 




3.5 


2.1 


3.7 








Xc0(l 3 P0) 


2.9 


1.39 


6.38 


2.5 


5.32 




3.1(8) 


2.90(59)(66)(30) 




X' C0 (2 3 P ) 


1.9 


1.11 
















Xc2(l 3 P 2 ) 


0.50 


0.44 


0.57 


0.28 


0.44 




0.45(8) 


0.31(5)(4) 


0.61(10) 


x' c2 (2 3 P 2 ) 


0.52 


0.48 
















%(!%) 


0.35 


0.22 


0.46 


0.46 


0.45 


0.17 










0.15 


0.11 




0.20 


0.21 


0.13 










0.10 


0.084 
















X60(1 3 P0) 


0.038 


0.024 


0.080 


0.043 












x' b0 (2 3 P ) 


0.029 


0.026 
















X62(1 3 P 2 ) 


0.008 


0.0056 


0.008 


0.0074 












X^(2 3 P 2 ) 


0.006 


0.0068 

















In Table 1 we compare our predictions for two-photon decay rates of heavy 
quarkonia calculated with the account of both relativistic and radiative cor rect ions 
with previous theoretical calculat ions and ex peri mental averages from PDcPl and 
very recent measurements by E83^^and Belief Collaborations. The application of 

rrni 

the quark-hadron duality to sum rules for two-photon decays yielded the following 
predictions T(r] c -> 77) S 7.4 keV, T( Xc q -> 77) < 6.1 keV, T( Xc 2 -> 77) < 3.9 keV. 

It follows from Eqs. (|14fl . (|15|1 that the ratio of two-photon decay rates of scalar 
and tensor states in the nonrelativistic approximation is independent of meson wave 
functions 

r NR ( 3 p -> 77) _ 15 

r NR(3p 2 _ 77 ) - 4 
Inclusion of radiative corrections l|17l) . i|18|) considerably increases this ratio 
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Indeed for the chosen values of a s we find i?^ R = 6.8 and i?S R — 5.5. On the other 

° cc bb 

hand, relativistic effects tend to decrease this ratio. The results of our calculations 
of the ratio Rqq = T( 3 P -> 7 7 )/r( 3 P 2 -> 77) using Eqs. (TTJ, CH, 03, CHI) are 
given in Table 2 in comparison with other theoretical predictions and experimental 
data. 
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Table 2. Ratio Rqq of two-photon decay rates of scalar ( 3 Po) and tensor ( 3 P2) 
states. Comparison of different theoretical predictions and experimental data. 





our 


Refill 


Ref|3 


R.eflTol 


ReflTTI 


ReflTTI 


CLEdiSI 


Rce(lP) 


5.8 


3.2 


11.2 


8.9 


12.1 


5.4(1.3) 


7.4(2.4)(5)(9) 


Rcc.{2P) 


3.7 


2.3 














4.8 


4.3 


10 


5.8 










4.6 


3.8 













Our analysis shows that both relativistic effects and radiative corrections play 
an important role in the investigation of two-photon annihilation rates of heavy 
quarkonia. Their inclusion considerably modifies theoretical predictions and brings 
them in accord with experimental data. For computations we use the relativistic 
quarkonium wave functions which were previously obtained in calculating the mass 
spectra of charmonium and bottomonium. The consistent account of the relativis- 
tic corrections to the decay amplitudes and the use of the relativistic quarkonium 
wave functions calculated in the framework of the relativistic quark model, which 
correctly predicts meson mass spectra and various decay rates, as well as inclusion 
of radiative corrections distinguish our approach from previous calculations. The 
increase of the experimental precision in measuring the charmonium two-photon 
decay rates and the observation of corresponding bottomonium decays will help to 
discriminate between different theoretical models. 
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